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ABSTRACT 



We investigate the ionized environment of the Seyfert 1 galaxy NGC 985 with a new 
Chandra-ttETGS observation and an archival BeppoSAX observation. Both spectra 
exhibit strong residuals to a single powerlaw model, indicating the presence of an 
Q^ , ionized absorber and a soft excess. A detailed model over the Chandra data shows that 

the 0.6-8 keV intrinsic continuum can be well represented by a powerlaw (T ~ 1.6) plus 
a blackbody component (kT=0.1 keV). Two absorption components are clearly required 
to fit the absorption features observed in the Chandra spectrum. The components have 
a difference of 29 in ionization parameter and 3 in column density. The presence of the 



i-^ ■ low ionization component is evidenced by an Fe M-shell unresolved transition array 

^ H (UTA) produced by charge states vn to xiii. The high ionization phase is required 

by the presence of broad absorption features arising from several blends of Fe L-shell 
transitions (Fe xvil-xxn). A third highly ionized component might also be present, 
but the data does not allow to constrain its properties. Though poorly constrained, 
the outflow velocities of the components (581 ± 206 km s _1 for the high ionization 
phase and 197 ± 184 km s _1 for the low ionization one) are consistent with each other, 
and with the outflow velocities of the absorption components observed in the UV. In 



addition, the low ionization component produces significant amounts of O vi, N V, and 
C iv which suggests that a single outflow produces the UV and X-ray features. The 
broad band (0.1-100 keV) continuum in the BeppoSAX data can be parameterized by 
a powerlaw (r w 1.4), a blackbody (kT=0.1 keV), and a high energy cutoff (E c ss 70 
keV). An X-ray luminosity variation by a factor 2.3 is observed between the BeppoSAX 
and Chandra observations (separated by almost 3 years). Variability in the opacity 
of the absorbers is detected in response to the continuum variation, but while the 
colder component is consistent with a simple picture of photoionization equilibrium, 
the ionization state of the hotter component seems to increase while the continuum 
flux drops. The most striking result in our analysis is that during both the Chandra 
and the BeppoSAX observations, the two absorbing components appear to have the 
same pressure. Thus, we suggest that the absorption arises from a multi-phase wind. 
Such scenario can explain the change in opacity of both absorption components during 
the observations, but requires that a third hotter component is pressure-confining the 
two phases. Hence, our analysis points to a 3-phase medium similar to the wind found 
in NGC 3783, and further suggests that such a wind might be a common characteristic 
in AGN. The pressure-confining scenario requires fragmentation of the confined phases 
into a large number of clouds. 



1. INTRODUCTION 

Highly ionized (or 'warm') absorbers 
(Halpern 1982) have been observed in 
about half of all UV and X-ray spectra of 
AGNs (Reynolds 1997, George et al. 1998, 
Crenshaw et al. 1999). The absorption 
lines of these systems are blueshifted with 
respect to the optical emission lines, which 
implies outflow of material. The mass loss 
rates inferred from these outflows can be 
a substantial fraction of (Mathur, Elvis & 
Wilkes 1995) or even larger than (Behar et 
al. 2003, Ogle et al. 2004) the accretion 
rates needed to power the AGN contin- 
uum, indicating that these 'quasar winds' 
play an important role in the mass and 
energy budgets of AGNs. 

The UV and X-ray absorbers must be 
closely related, as exactly the same AGNs 
show both (Mathur, Elvis, & Wilkes 1995; 
Mathur, Wilkes, & Elvis 1998; Crenshaw 
et al. 1999; Monier et al. 2001). Further- 
more, their kinematic properties also ap- 
pear to be related, as both absorbers show 
similar outflow velocities (~1000 km s _1 ). 
Given the intrinsic complexity of UV and 
X-ray spectra, along with the large differ- 
ence in spectral resolution in the two bands 
(with an order of magnitude lower spec- 
tral resolution in the X-ray region), it has 
been extremely difficult to connect these 
absorbers in detail, leading to highly com- 
plicated models to describe them. Never- 
theless, recent studies using more compre- 
hensive models have shown that the obser- 
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vational data are certainly pointing to a 
common nature for the UV and X-ray out- 
flows (see §5.1). 

New studies based on high quality data 
have shown that the absorbers can be de- 
scribed in a surprisingly simple way. For 
example, Krongold et al. (2003) mod- 
eled the 900 ksec Chandra HETG expo- 
sure on NGC 3783 (Kaspi et al. 2002) 
and found that most of the absorption fea- 
tures present in the spectrum could be re- 
produced by a model with only two ab- 
sorption components, and that these com- 
ponents appeared to be in pressure bal- 
ance. This result was confirmed by Net- 
zer et al. (2003), who found a third highly 
ionized component, also in pressure equi- 
librium with the other two. Is this single 
object result only a coincidence? Or, on 
the contrary, does it suggest that promis- 
ing advances in our understanding of the 
physical conditions of AGN winds can be 
achieved, carrying important repercussions 
in our understanding of the structure of 
quasars (e.g. Elvis 2000)? 

Here, we present Chandra HETG and 
BeppoSAX observations of the Seyfert 1 
galaxy NGC 985 and show that the X-ray 
spectra of this source can also be well de- 
scribed by a simple model with two absorp- 
tion components. Furthermore, as for the 
ionized absorber toward NGC 3783, these 
components are consistent with pressure 
equilibrium, with one phase having suffi- 
ciently low ionization to produce substan- 
tial absorption in the UV region. We fur- 
ther discuss the implications of these re- 
sults. 
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1.1. NGC 985 

NGC 985 (Mrk 1048) is a Seyfert 1 
galaxy located at redshift 0.04274±0.00005 
(12814±15 km s -1 ; Arribas et al. 1999, 
based on stellar absorption features). NGC 
985 is a peculiar galaxy with a promi- 
nent ring shaped zone at several kilopar- 
secs from the nucleus (de Vaucouleurs & 
de Vaucouleurs 1975), suggesting that the 
galaxy is undergoing a merging process. 
Further evidence for a merger has been 
found through optical and near-IR ob- 
servations, which have revealed the pres- 
ence of a double nucleus (Perez Garcia & 
Rodriguez Espinosa 1996, and references 
therein). The nuclei are separated by 3" 
(2.5 kpc of projected linear distance; Ap- 
pleton et al. 2002), which indicates that 
the intruder galaxy responsible for the for- 
mation of the ring is sinking in the nuclear 
potential of the primary galaxy (Perez- 
Garcia & Rodriguez-Espinosa 1996). As 
in other galaxies undergoing a merger, 
the IR luminosity of NGC 985 is extreme, 
Lm = 1.8 x 10 11 Lq, which puts this object 
in the LIRG group. 

An ionized absorber has been observed 
in the X-ray spectra of NGC 985. This 
source has been observed in the X-rays 
with different satellites and is a strong X- 
ray source. The first detection was with 
the Einstein Observatory, measuring an 
0.2-4.0 keV flux of 1.76xl0 _11 erg s" 1 cm" 2 
(a luminosity of 1.4 xlO 44 erg s -1 ; Kruper, 
Urry, & Canizares 1990). Subsequently, 
spectral analysis of ROSAT-PSPC data 
revealed complex features (Brandt et al. 
1994). A simple power law could not fit 
the PSPC spectrum of NGC 985, and the 
presence of a warm absorber, as well 
continuum soft excess, was suggested. In 



1996, NGC 985 was observed by ASCA 
with a total exposure of 40 ksec. The anal- 
ysis of these data confirmed the presence of 
an ionized absorber in the spectra of this 
source (Nicastro et al. 1998, 1999). Nicas- 
tro et al. found that photoelectric bound- 
free edge absorption could not account for 
all the features observed, and suggested 
that resonant bound-bound line absorption 
played an important role. 

This possibility was supported by HST- 
STIS and FUSE ultraviolet observations of 
the nucleus in NGC 985, which have re- 
vealed the presence of six narrow absorp- 
tion components detected in Ly a, C IV, 
N v and O vi (Arav 2002). Five of these 
systems are observed in outflow, with ve- 
locities (with respect to the centroid of the 
narrow emission lines) of 780, 670, 519, 
404, and 243, km s _1 (hereafter compo- 
nents 1 to 5). The sixth component is 
observed inflowing at —140 km s~ x . As 
demonstrated by Arav, small to moderate 
saturation is present in the Ly a trough 
of these components. Once the effects of 
saturation are accounted for, at least com- 
ponent 4 is consistent with gas that also 
produces X-ray absorption, i.e. a highly 
ionized or 'warm' absorber, with a column 
density of 10 21,4 cm -2 . 

In this paper we present 2 moderate ex- 
posure observations of NGC 985 (~ 100 
ksec each) with Chandra and BeppoSAX. 
In §2 we describe the reduction and spec- 
tral analysis of the Chandra-HETGS data. 
The analysis was carried out in Sherpa 1 
(Freeman, Doe, & Siemiginowska 2001), 
with the inclusion of the code PHASE 
(Krongold et al. 2003) to model in a 
self-consistent way the ionized absorber 
present in the spectra. In §3 we present the 
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BeppoSAX data, and in §4 we analyze the 
variability effects of the source on the ab- 
sorber. In §§5 and 6 we discuss the differ- 
ent scenarios consistent with the absorber 
and their implications. 

2. THE CTmndra-HETG SPECTRUM 

On 2002 June 24, the High Energy 
Transmission Grating Spectrometer (HETGS, 
Canizares et al. 2000) on board the 
Chandra X-ray Observatory (Weisskopf et 
al. 2000) observed NGC 985 (see Table 
1). The observation was carried out with 
the Advanced CCD Imaging Spectrometer 
(ACIS; Garmire et al. 2003) in the focal 
plane, and had a duration of 80 ksec. We 
processed the event files of this observa- 
tion with the Chandra Interactive Analysis 
of Observations (CIAO 2 ; Fruscione 2002) 
software (Version 2.3), following the on- 
line data analysis "threads" provided by 
the Chandra X-ray Center (CXC 3 ), and 
extracted source and background 1st or- 
der spectra (positive and negative orders) 
and responses for the High-Energy Grat- 
ing (HEG) and Medium-Energy Grating 
(MEG) configurations. We accounted for 
the ACIS efficiency degradation (Marshall 
et al. 2003) using the script contamabs 4 
(version 1.1). The net exposure of the ob- 
servation is 77.7 ksec and the total num- 
ber of counts in the dispersed spectrum is 
13,979. 

We then combined the -1st and +lst or- 

1 http: / / cxc.harvard.edu/sherpa 
2 http://asc. harvard.edu/ciao 

3 http://asc. harvard.edu/ciao/documents_threads. html 

4 http://space. mit.edu/CXC/analysis/ACIS_ Con- 
tam/ACIS_ Contam.html 



der spectra to improve the signal to noise 
(S/N) ratio of the data. Similarly, after 
checking for wavelength scale consistency, 
the HEG and MEG spectra and their re- 
sponses were coadded in the overlapping 
region, to further increase the S/N ratio. 
The S/N ratio per MEG resolution element 
(«0.02 A) varies from 7 to ~ 2 in the 
2-20 A range, a factor « 6 lower than the 
highest S/N ratio HETG spectrum of a Sy 
1 galaxy (NGC 3783; Kaspi et al. 2002). 
All the spectra presented throughout the 
paper (except where otherwise stated) cor- 
respond to a binning factor of 20 (i.e we 
binned the spectrum in bins of size 0.1 
A, about 5 times the MEG resolution el- 
ement), giving a S/N ratio ranging from 
14 to 3 per bin (in the 2-20 A range). The 
spectrum is shown in Figure la. 

2.1. Spectral Analysis 

In the calculations reported here, we at- 
tenuated the continuum by an equivalent 
hydrogen column density of 3.0 x 10 20 cm -2 
to account for the Galactic absorption by 
cold gas in the direction to the source 
(Starck et al. 1992). We have explored 
only photoionization equilibrium models, 
and we have assumed solar elemental abun- 
dances (Grevesse & Noels 1993). 

2.1.1. Fitting the Continuum 

To model the intrinsic continuum of the 
source, we first used a simple power law. 
This model could not fit the data over the 
entire energy range (0.6-8 keV, 1.5-20 A), 
as can be seen from the residuals (Fig. 2a). 
The model shows both positive and neg- 
ative deviations (with significance > 2a) 
typical of the ionized absorber and soft ex- 
cess commonly observed in Seyfert 1 galax- 
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ies (e.g. Reynolds 1997). 

We repeated the power law fit using 
only data between 3 and 8 keV (we also 
excluded the Fe Ka line region around 
6.4 keV). At these energies, absorption by 
moderately large columns of ionized gas 
(i.e. < 10 23 cm -2 ) or low columns of neu- 
tral gas (i.e. < 10 21 cm -2 , Galactic gas) 
do not modify significantly the shape of 
the intrinsic continuum (e.g. Krongold et 
al. 2003). The results of this model are 
summarized in Table 2, and the residuals 
(extrapolated to the whole spectral range) 
are presented in Figure 2b. The presence 
of strong positive and negative deviations 
below 3 keV in this model shows that both 
an ionized absorber and a soft excess are 
likely to be present. To account for the 
presence of the excess, we further included 
a blackbody component in our continuum 
determinations, and the positive residuals 
decreased significantly (see Fig. 2c). How- 
ever, the model still shows strong nega- 
tive residuals at energies characteristic of 
an ionized absorber. 

Finally, we fitted the continuum and 
the absorption together (see §2.1.2 for the 
details of the ionized absorber). In this 
model we left all the continuum compo- 
nents (power law and blackbody) free to 
vary, and we assumed that the absorbing 
gas also covers the soft excess. The re- 
sults obtained for the power law with this 
approach are in excellent agreement with 
the results obtained when considering only 
the 3-8 keV range. Table 2 presents the 
values for the best continuum fit. Figure 
la shows the fluxed data over the entire 
spectral range (1.5-20 A, 0.6-8 keV), to- 
gether with (a) the intrinsic continuum of 
the source, (b) the continuum attenuated 



by Galactic absorption, and (c) the ob- 
served continuum (i.e. the continuum fur- 
ther attenuated by the bound-free transi- 
tions from the ionized absorber, see details 
in Fig. caption). 

2.1.2. Modeling the Ionized Absorber 

We ran PHASE with three free parame- 
ters for each component necessary to fit the 
ionized absorber: (1) the ionization param- 
eter (defined over the entire Lyman contin- 
uum, as U — Q(H)/4irr 2 n H c with Q(H) 
being the rate of H ionizing photons, r the 
distance to the source, n# the H density, 
and c the speed of light.), (2) the equiva- 
lent hydrogen column density N#, and (3) 
the outflow velocity. Other parameters of 
the code that we kept fixed are the internal 
micro-turbulent velocity of the absorbers 
and the intrinsic Spectral Energy Distri- 
bution (SED) of the source (Fig. 3a). The 
turbulent velocity of the gas is very diffi- 
cult to constrain, since (1) single absorp- 
tion lines in the spectrum are unresolved 
by the HETGS (which has a FWHM res- 
olution of ss 400 km s _1 at 17 A) and (2) 
most of the observed features are blends 
of several transitions. We set this parame- 
ter to 350 km s _1 , the value most suitable 
to the data after several tests comparing 
the measured equivalent widths (EWs) of 
a few unblended lines (Si xiv, Si xill, Mg 
XII, Mg XI, Fe xvn) with the EWs pre- 
dicted by our model. 

To approximate the SED, we proceeded 
as follows: in the X-ray region we used 
the power law T, as well as the blackbody 
contribution, as inferred from the fits to 
the continuum presented in §2.1.1. To ob- 
tain the exact parameters for the black- 
body component several iterations were 
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performed. The BeppoSAX data of NGC 
985 show the presence of a high energy cut- 
off in the intrinsic continuum (see §3.1.1), 
which we modeled with a break in the X- 
rays at around 70 keV, and a powerlaw 
with T = 5 above this energy. Based 
on HST-STIS observations of NGC 985, 
Bowen, Pettini & Blades (2002) found a 
flux of 2.54 xl0~ 26 erg s _1 cm -2 hz -1 at 
1230 A. In order to determine the flux at 
the Lyman limit, we extrapolated to 912 
A assuming a power law with photon in- 
dex T =2. The flux at 912 A turned out 
to be 1.83 xl0~ 26 erg s _1 cm -2 hz -1 . Ex- 
trapolation of the X-ray continuum down 
to the Lyman limit gives a flux two or- 
ders of magnitude smaller than measured 
by HST, which indicates the presence of a 
soft X-ray-UV energy break, with a steep- 
ening of the spectrum between the two 
bands (Fig. 3a). We assumed that this 
break occurs at 0.1 keV. The actual en- 
ergy of this break is unknown (due to the 
unobservable region in the UV- X-ray range 
associated with Galactic absorption). We 
restricted the SED at low energies using 
(not simultaneous) data from NED 5 (see 
Fig. 3a). For simplicity we extended the ex- 
trapolation between 1230 A and 912 A up 
to 100 /im (i.e. we considered a power law 
with T =2 from 912 A to 100 /im). Below 
100 /im we introduced a low energy cut- 
off with T = —3.5. The final SED used 
to model the Chandra data is presented in 
Figure 3a. The differences in the infrared 
region between NED data and our assumed 
SED do not have any effect on our analysis 
(the only effect such differences produce is 
an increase of the Compton temperature 
by a factor <1.1, when the NED data are 



considered) . 

2.1.3. The Ionized Absorber 

As is the case for other ionized absorbers 
in AGNs (Kaspi et al. 2001; Blustin et. al 
2002; Krongold et al. 2003; Netzer et al. 
2003 for NGC 3783; Kaastra et al. 2002 
for NGC 5548; Steenbrugge et al. 2003 for 
NGC 4593; Blustin et al. 2003 for NGC 
7469) the ionized absorber in the spectrum 
of NGC 985 could not be well fitted by 
a single absorbing component (x 2 /d.o.f.= 
201.6/177). Therefore, we attempted to 
fit the spectrum with the next simplest 
model, by including a second component. 
In this case, we obtained an acceptable fit 
to the data (x 2 /d.o.f.= 179.03/174; an F- 
test indicates that the second component is 
required at a significance level larger than 
99.9%). Table 3 lists our results, and the 
best fitting model is presented in Figure 4 
(and also in Fig. lb, where we show the 
model for the full spectral range plotted 
over the fluxed spectrum) . Most of the nar- 
row features observed in Figure 4 have a 
significance < 2cr, but still the absorber is 
clearly detected and constrained through 
the presence of relatively broad features 
due to both photoelectric bound-free ab- 
sorption and blends of resonant bound- 
bound transitions (i.e. Fe L-shell and Fe 
M-shell absorption). To highlight the pres- 
ence of these features (which are significant 
at a level > 2.5<r), in Figure 5 we present 
our model and the data in the 12-18 A 
range, with a binning; size of 0.15 A (b in- 
ning factor of 30). 

Two different ionization absorbing com- 

5 NASA/IPAC Extragalactic Database, 
http: / / nedwww.ipac.caltech.edu 
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ponents are clearly required by the data. 
The presence of significant features pro- 
duced by both Fe L-shell transitions (Fe 
xvil-xxil) and Fe inner M-shell transitions 
(Fe vii-xiii, the unresolved transition ar- 
ray or UTA) is incompatible with a sin- 
gle ionization degree for the absorbing gas. 
This is shown in Figure 6, where the con- 
tribution from the two components to the 
absorption is presented separately. The 
high ionization phase (hereafter HIP), with 
logU=1.34, gives rise to the absorption by 
Oviii, Ne ix-x, Mg xi-xii, Si xm-xiv, and 
Fe L-shell. The low ionization phase (here- 
after LIP), with logU=-0.12, is responsible 
for the absorption by charge states O VII, 
Ne v-ix, Mg IX, and the Fe M-shell UTA. 
The components are well separated in ion- 
ization degree, the HIP has a value « 29 
times larger in U than the LIP. Accord- 
ingly, the HIP has a higher temperature 
(logT=5.8 [K]) than the LIP (logT=4.5 
[K]) by a factor of 21. These temperatures 
were calculated assuming photoionization 
equilibrium in each component. Figure 7 
shows the theoretical spectral features pre- 
dicted for the absorption components. Ta- 
ble 5 presents the predicted column densi- 
ties for the dominant charge states in the 
LIP and the HIP. 

As has been observed in other ionized 
absorbers (see above), the HIP has a larger 
equivalent H column density than the LIP, 
in the case of NGC 985 by a factor of 
3. The outflow velocity of the compo- 
nents is 581±206 km s" 1 for the HIP and 
197±184 km s" 1 for the LIP. These veloc- 
ities, though ill-determined, are consistent 
with each other (within ~ la) and with 
those observed in the HST and FUSE sys- 
tems (see Fig. 8 for more details). A fit 



fixing the velocity of both HIP and LIP 
at the velocity of the strongest UV com- 
ponent (4, outflowing at 404 km s -1 ; see 
§1 and Fig. 8) yields an acceptable fit to 
the data, with parameters fully consistent 
with those obtained before. 

In the former analysis we have assumed 
that the soft excess detected in the data 
is covered by the absorber. Soft X-ray ex- 
cesses have been studied in low resolution 
spectra for more than a decade (e.g., Elvis 
et al. 1991, and references therein) and 
they have been found to vary on timescales 
of tens of days. This effect has been con- 
firmed in high resolution spectra (Netzer 
et al. 2003), and suggests that the bulk 
of the soft X-ray emission comes from a 
compact region (less than 10 16 cm across; 
Elvis et al. 1991). It has been suggested 
that the emission might have a thermal na- 
ture and might originate in the inner edges 
of accretion disks (e.g. Czerny & Elvis 
1987). Therefore, it is reasonable to as- 
sume that the absorbing material lies fur- 
ther away from the central engine than the 
soft X-ray emission. Nevertheless, we note 
that the data does not allow to distinguish 
whether the soft component is covered by 
the absorber or not. A fit assuming an un- 
covered excess is statistically indistinguish- 
able from a fit which assumes full cover- 
ing. If the excess is not absorbed, then 
the column density of the LIP is 1.6 times 
larger and the normalization of the black- 
body 2.3 times smaller. The rest of the 
parameters in this model are fully consis- 
tent with the values obtained before. We 
stress that whether the soft excess is cov- 
ered or not covered by the absorber, this 
has a negligible effect on the conclusions 
presented in this paper. 
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We find no significant evidence of an 
ionized emitter in the spectrum. Resonant 
and forbidden Ovil lines have a statistical 
significance < 1 sigma. We estimate upper 
limits for the equivalent width of 218 mA 
for the resonant line and 185 mA for the 
forbidden line. Although an emitter might 
be present, better data are required for its 
detection. 

2.2. Constraining the Ionized Ab- 
sorber 

Despite the limited S/N ratio of the 
Chandra spectra, the physical properties 
of the ionized absorber in NGC 985 can 
be constrained. As demonstrated by Be- 
har, Sako, & Kahn (2001) and Krongold 
et al. (2003), the UTA is a robust indica- 
tor of the ionization level of the gas, as the 
general shape and position of this feature 
strongly depend on U for a given ioniza- 
tion component. The 12-15 A Fe L-shell 
absorption complex also helps to limit the 
physical state of the gas. 

To illustrate this point, we have con- 
structed a series of models with ionization 
parameters different from those deduced 
for our best fitting model (see §2.1.3 and 
Table 3). The comparison of these models 
with the data reveals clear discrepancies 
(see Fig. 9). The first model we considered 
is presented in Figure 9a, and consists of 
a factor 3 less ionized HIP (logU=0.8). In 
this model, lower charge states of Fe would 
be significantly present in the HIP (e.g. Fe 
xv-xvi), and these would produce abun- 
dant Fe M-shell UTA absorption at wave- 
lengths between 15 and 16 A (see Fig. 9a). 
This is clearly not observed in the data. 
At the same time, the deep absorption fea- 
ture observed in the data around 12.85 A 



(arising from Fe XIX-XX), is missing in the 
model. Thus, the HIP cannot contain less 
ionized gas. 

The ionization state of the LIP is also 
limited by the data, particularly by the 
UTA. Figure 9b shows the comparison of 
the data with a second model that includes 
a factor 4 more ionized LIP (logU=0.5). 
This model shows an excess of absorp- 
tion at 15- 16 A produced by Fe xiv- 
XV, and a lack of absorption at longer 
wavelengths (due to a diminished fraction 
of Fe in charge states vn-x). A third 
model consisting of a factor 5 less ionized 
LIP (logU=-0.8) is presented in Figure 9c. 
This model is inconsistent with the data 
also. A low value of U shifts the UTA to 
higher wavelengths (peak of absorption at 
17 A) because an important fraction of Fe 
is now present in charge states iv-vi. At 
the same time, the lack of Fe x-xil pro- 
duces an underestimation of the absorption 
around 15.9-16.4 A. From this analysis is 
evident that even lower values of U for the 
LIP (logU<-0.8) would produce even more 
shifted UTAs (produced by even lower Fe 
charge states), and therefore, higher dis- 
crepancies between data and model. Thus, 
the data impose important restrictions on 
the ionization state of the absorbing gas. 

A third absorber with an intermediate 
value of the ionization parameter (between 
the HIP and the LIP) and high H col- 
umn density can also be ruled out. Fig- 
ure 9d illustrates the inconsistencies that 
arise when a third absorber with logU=0.6 
and logNtf=21.5 [cm -2 ] is included in our 
model. Such an absorber severely overpre- 
dicts the absorption by the UTA from 15 
to 16 A (because of a high contribution 
to the absorption in the model by Fe xiv- 
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xvi) . However, it is not possible to rule out 
the presence of additional absorbing com- 
ponents with lower H column density. 

To quantify the limits imposed by the 
data to additional absorbing material we 
included a third absorption component in 
our model and refit the data. In Figure 10a 
we show the logU-logN^ confidence region 
for this hypothesized absorber. The plot 
shows that gas with U between the HIP 
and the LIP and logN^ > 21.1 [cm" 2 ] (half 
the column density of the LIP and 20% 
that of the HIP) can be ruled out (at the 
3<T level). Figure 10a also shows that the 
presence of material less ionized than the 
LIP and logN# > 21 [cm -2 ] can also be 
ruled out (at 3a significance). 

On the other hand, a third component 
with very high ionization parameter and 
considerable column density can be present 
in the spectra of NGC 985. Figure 10b 
presents the confidence region for a 3rd 
component consisting of high column den- 
sity and high ionization gas. As observed 
in the figure such a component is not well 
constrained by the data, but in order to 
have larger column density than the HIP 
(logNff >21.8 [cm" 2 ]) requires logU>2 (i.e 
5 times larger U than the HIP). In fact the 
discrepancies (with significance < 2a) be- 
tween 1.6 and 1.7 A, and near 10.6 A (see 
Fig. 11) hint at the presence of such mate- 
rial, as they are consistent with absorption 
by Fe XXIV- XXVI. 

It is clear, then, that the gas tends to 
clump in two different degrees of ioniza- 
tion, with the majority of the low ionized 
material around the LIP. A third, very high 
ionization component is hinted at by the 
data, but better S/N ratio data are needed 
to establish its existence and constrain its 



physical properties. 

2.2.1. The Effects Induced by the Lack of 
Accurate Dielectronic Recombina- 
tion Rates 

Our analysis of the LIP is based mostly 
on the Fe M-shell UTA, and low tem- 
perature Dielectronic Recombination Rate 
(DRR) coefficients have not been calcu- 
lated or measured yet for these ions. Kron- 
gold et al. (2003) and Netzer et al. (2003) 
reported that the models were unable to 
simultaneously reproduce the absorption 
by the UTA and two lines of Si (Six at 
A6.850 and Sixi at A6.775). Netzer et al. 
(2003) proposed that this was an effect in- 
troduced by the lack of these DDR. This 
means that, for a given value of the ion- 
ization parameter, there is a systematic 
shift in the ionization balance of Fe (to- 
ward higher charge states) with respect to 
the rest of the species. Netzer (2004) and 
Kraemer, Ferland, & Gabel (2004) have 
estimated these rates. Their results show 
that there is indeed an increase in the ion- 
ization parameter at which the Fe M-shell 
ions peak, offering a solution for the mis- 
match between the UTA and Si lines in the 
spectrum of NGC 3783. 

In the spectrum of NGC 985, besides 
the UTA, there are no other significant ab- 
sorption features that allow us to constrain 
the ionization parameter of the LIP. Ac- 
cepted at face value, this means that we are 
underestimating the ionization parameter 
(and temperature) of the LIP, because of 
the lack of Fe M-shell DRR. However, from 
our previous analysis on NGC 3783 (Kro- 
ngold et al. 2003), we estimate that the 
error cannot exceed a factor of 2 in U (i.e. 
Alogf/ < 0.3), a factor consistent with 
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the one obtained by Netzer et al. (2003) 
and Kramer, Ferland, & Gabel (2004). We 
have repeated our analysis shifting up by 
a factor of 2 the ionization parameter of 
the LIP in all elements except Fe. With 
the quality of the present data, this dis- 
tinction does not modify our results in any 
significant way. In fact, estimates with this 
new approach of the column density and 
ionization parameter of possible additional 
absorbing components yield to limits fully 
consistent with those presented in the for- 
mer section, showing again that the gas 
clumps in two, or maybe three, absorbing 
components. 

While our main conclusions will not 
change due to this underestimation of U, 
this effect should be kept in mind. There- 
fore, we have set the positive error bar in 
the ionization parameter of the LIP to be 
consistent with a factor of 2 change in U, 
in order to reflect the underestimation in- 
troduced by the lack of low temperature 
DRR for Fe. 

Recent calculations by Gu et al. (2003), 
and measurements by Savin et al. (2002a, 
2002b, 2003) have produced more accu- 
rate values for Fe L-shell DRR. At tem- 
peratures close to the one derived for the 
HIP ~ 10 6 K (~ 86 eV) the differences be- 
tween the old estimates and the new ones 
are relatively small (although these differ- 
ences become larger at lower temperatures, 
i.e. 10 5 K). Thus, although the new calcu- 
lations have not been incorporated into the 
photoionization codes, we expect the un- 
certainties introduced into the HIP model 
by the uncertainties of the Fe L-shell DRR 
currently used in the codes to be within 
the error bars derived in our analysis. In 
fact, our model for this component is able 



to fit the Fe L-shell lines, along with K- 
shell lines arising from other elements, like 
Si, Mg, and Ne (see Fig. 4). Although the 
significance of these narrow lines is limited, 
it is clear that the ionization state of the 
gas is consistent with the observed charge 
states. We note that the same was found 
for the high S/N ratio spectrum of NGC 
3783 (Krongold et al. 2003; Netzer et al. 
2003). Therefore, we do not expect a sys- 
tematic underestimation of U in the HIP. 

3. BeppoSAX DATA 

To further study the nature of the ion- 
ized absorber, in particular to constrain 
the absorber response to the changes in 
the continuum of the central source, we 
retrieved and reprocessed the previously 
unpublished NGC 985 data obtained by 
BeppoSAX (Boella et al. 1997; Piro et 
al. 1997) on 1999 August 29 (see Table 
1). The data have a total MECS expo- 
sure of 90 ks. We followed the standard 
reduction procedure (Fiore et al. 1999) for 
the LECS, MECS, and PDS data. We ex- 
tracted LECS and MECS source spectra 
from circles of radius 8' and 4' respectively, 
and the same extraction regions were used 
to extract background spectra from the 
LECS and MECS "blank fields." The PDS 
instrument has no imaging capability and 
is a collimated detector with a field of view 
(FOV) of 1°.3 (FWHM). It makes use of 
rocking collimators for background moni- 
toring. Using the ROSAT All Sky Survey, 
we searched the FOV of the PDS in the 
direction to NGC 985, and did not find 
any other source in this region capable of a 
significant contamination of the PDS spec- 
trum (Mrk 1044 is the closest bright X-ray 
source to NGC 985 in the sky, at an angu- 
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lar distance of 1.2°, but it was outside the 
FOV of the PDS during the observation). 

Modest variability was observed dur- 
ing the BeppoSAX observation of NGC 
985. In the 0.1-2.0 keV band (LECS data), 
source variations were detected up to a fac- 
tor 1.45 on a time-scale of about 50 ksec. 
Variability in the 2.0-10 keV range (MECS 
data) was detected up to a factor of 1.3 in 
the same time-scale. The quality of the 
data does not allow for a separation in low 
and high state spectra, therefore we mod- 
eled the full time integrated dataset. 

3.1. Spectral Analysis 

3.1.1. Fitting the Continuum 

To model the broadband (0.1-100 keV) 
BeppoSAX continuum of NGC 985, we 
first followed the same approach used with 
the Chandra data, constraining the intrin- 
sic powerlaw in the 3-8 keV range. The 
residuals of the extrapolation of this model 
to the whole BeppoSAX band are pre- 
sented in Figure 12. The presence of both 
an ionized absorber and a soft excess is ev- 
ident, together with that of a high energy 
cutoff at energies larger than 60 keV. We 
then refitted the data including an expo- 
nential cutoff to model the continuum at 
high energies, as well as a blackbody and 
an ionized absorber (see §3.1.2) to model 
the soft excess and the absorption at low 
energies. 6 Results of our best continuum 
fits are shown in Table 2. The intrinsic 
continuum of the source further attenuated 
by Galactic absorption is shown in Figure 
13a (dashed line). 

6 The inclusion of a cold reflector, often observed 
in the high energy spectra of Seyfert galaxies (e.g. 
NGC 3783, Kaspi et al. 2001, de Rosa et al. 2002), 



The BeppoSAX data indicate a flatter 
power law photon index (Ar ~ 0.2 ±0.07) 
and a larger flux in both the power law 
and the blackbody components (by a factor 
w 2.3) with respect to the Chandra data. 
These differences are reflected in the SED 
used to fit the ionized absorber in the Bep- 
poSAX data (see Fig. 14). 

3.1.2. Modeling the Ionized Absorber Present 
in the BeppoSAX Data 

Modeling of absorption features in low 
resolution data is extremely difficult be- 
cause the lack of detail does not allow 
the determination of the exact ionization 
state of the gas and the number of ab- 
sorbing components. Because of this, we 
decided to use the solution found for the 
high resolution Chandra data, to constrain 
the physical properties of the gas at low 
resolution. We assumed that two compo- 
nents contribute to the absorption, as in- 
ferred from the Chandra data, and fixed 
the outflow velocity and the turbulent ve- 
locity of the absorbers to the values de- 
duced in §2.1.3 (these parameters are es- 
sentially unconstrained when modeling low 
resolution data). 

First, we left the ionization parameter 
and the column density of both compo- 
nents free to vary while fitting the data. 
An acceptable fit was obtained (reduced 
X 2 = 1.063 for 112 degrees of freedom), 
with column density values in both absorp- 
tion components closely consistent with 
the values obtained modeling the Chan- 
dra data: we find column density values of 

is not required by the BeppoSAX data of NGC 
985. We find a 2a upper limit for the amount of 
reflection R< 0.7, see Table 2. 
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logN#=21.95±.24 [cm" 2 ] (AlogN ff =0.14 
[cm" 2 ]) for the HIP and logN# =21.38±.26 
[cm" 2 ] (AlogN H =0.01 [cm" 2 ]) for the LIP. 
Hence, we next froze the column densities 
to the values obtained with Chandra, and 
fit the data letting only the ionization pa- 
rameters vary. The results are reported in 
Table 4, and are discussed in the next sec- 
tion. Figure 13 shows the data and the 
best fit model and Table 5 shows the pre- 
dicted column densities for the dominant 
charge states. 

4. VARIABILITY ANALYSIS 

4.1. Intensity and Spectral Changes 
in the Continuum 

A comparison of the 2002 Chandra and 
1999 BeppoSAX observations indicates a 
flux drop by a factor of 2.3 in the 0.1-10.0 
keV range. Variability is observed in both 
the soft X-ray (0.1-2.0 keV) and hard X- 
ray (2.0-10.0 keV) bands. Since the hard 
X-ray region is barely affected by the ion- 
ized absorber, the change of flux is indica- 
tive of intrinsic source variations. Our con- 
tinua fits indicate a steepening in slope of 
Ar « 0.2 ± 0.07. Any attempt to fit both 
datasets with a single photon index value 
yields unsatisfactory results. For instance, 
we find a x 2 = 59.1 for 50 degrees of free- 
dom when fitting (between 3 and 8 keV) 
the Chandra data with a value of T frozen 
to 1.4, this imply a A% 2 = 12.4 with re- 
spect to the best fitting model (an F-test 
gives a significance larger than 99.9% for 
this difference). Fitting the BeppoSAX 
data with a photon index equal to 1.6 pro- 
duces even worse results. This difference 
could be indicating a real slope variation 
or could be the result of cross-calibration 



effects between the two instruments. We 
will follow the analysis assuming T=1.4 for 
the BeppoSAX data, noticing that this dif- 
ference does not have a significant effect on 
our conclusions. 

To investigate the significance of the 
variability of the soft emission component, 
we attempted to fit the BeppoSAX data 
while freezing the blackbody component 
used to parameterize the soft excess to the 
values obtained for our best fitting Chan- 
dra model. Every model tried imposing 
this constraint led to an underestimate of 
the continuum level below 1 keV, pointing 
to a genuine variation of the flux of the soft 
component between the two observations. 
We note from Table 2 that only a change in 
the intensity of the emission is required by 
the data, as the best fit solutions to both 
datasets are consistent with a blackbody 
temperature of 0.1 keV. 

4.2. Changes in the Opacity of the 
Ionized Absorber 

The ionization degree of the ionized ab- 
sorber along the line of sight to NGC 
985 has also changed between the Chan- 
dra and BeppoSAX observations. An F- 
test comparing the opacities deduced from 
the Chandra best fitting model with those 
for the BeppoSAX best fitting model shows 
that the difference is statistically signifi- 
cant at the 99% confidence level. To better 
quantify the changes in the ionized mate- 
rial, it is convenient to introduce a different 
ionization parameter, U x , defined through 
the number of ionizing photons between 
0.1 and 10.0 keV (Netzer 1996). This def- 
inition makes the changes in the ioniza- 
tion parameter (now relevant only to the 
X-ray species) directly proportional to the 
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change in the X-ray flux, as opposed to 
bolometric flux changes, which are not di- 
rectly observable. Thus, it is preferable in 
studies of the effects of flux variability over 
ionized gas in the X-ray regime. The two 
different definitions of the ionization pa- 
rameter scale as logU :E =logU - 2.43 for the 
Chandra SED, and logU x =logU - 2.08 for 
the BeppoSAX SED. 

The best fit to the Chandra data gives 
a value of logU x =-2.55±0.09 for the LIP. 
If the gas is in photoionization equilibrium 
and is not physically confined, a change of 
flux by a factor 2.3 would induce a change 
in the X-ray ionization parameter by the 
same factor, yielding a value of logU x =- 
2.20. The best fit value of logTX, for the 
BeppoSAX data is -2.30±0.10, consistent 
with the predicted value. We conclude that 
the LIP is responding to the changes in 
the continuum and thus is probably both 
photoionized and in equilibrium. 

The value of logU x for the HIP in the 
Chandra data model is -1.09±0.10. If the 
absorber is in equilibrium, the flux change 
observed between observations implies a 
value of logU a; =-.73 for the HIP present in 
the BeppoSAX spectrum. The value ob- 
tained in our best model for the BeppoSAX 
data is logUa;=-1.25±0.13. These two val- 
ues of the ionization parameter are clearly 
inconsistent with each other. Thus, this 
component has indeed changed ionization 
state between the two observations, but 
does not seem to be following the changes 
on the continuum. On the contrary, the 
degree of ionization of this gas seems to 
have increased while the flux decreased. 

Such behavior has several possible ex- 
planations. (1) The gas is out of equi- 
librium: the photoionization equilibrium 



time scale is inversely proportional to the 
density in the gas during both increasing 
and decreasing flux intensity phases (e.g. 
Nicastro et al. 1999) . If the density is 
low enough, this timescale can be larger 
than the variability timescale of the central 
source, and the medium does not have time 
to adapt to the changes in the continuum. 
In such conditions, delays in gas responses 
may lead to scenarios where the ionization 
degree of the gas decreases while the flux 
increases. Such cases have been already 
found in other ionized absorbers around 
AGN, for instance NGC 4051 (Nicastro et 
al. 1999). (2) The gas is not in photoion- 
ization equilibrium because another mech- 
anism, collisional ionization (due, for in- 
stance, to shock heated gas) dominates. 
In this case, obviously, changes in the gas 
opacity are essentially unrelated to changes 
in the source flux. (3) The gas is indeed in 
photoionization equilibrium, and has un- 
dergone a significant change in its density. 
This is possible if the gas is pressure con- 
fined by an external medium. (4) Finally, 
it is possible that the absorbing gas is not 
the same during the two observations. X- 
ray absorbers can change independently of 
flux variations (see George et al. 1998) be- 
cause of transverse motion of the material 
or because of condensation and dissipation 
of clouds out of a hot medium (Krolik & 
Kriss 2001). In the following sections we 
will explore further the validity of these 
different scenarios. 

5. ABSORBER PHYSICS DISCUS- 
SION 

The ionized absorber present in NGC 
985 is suggestive of a simple picture, being 
constrained by two strong absorption com- 
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ponents. The same situation was found for 
the ionized absorbers in IRAS 13349+2438 
(Sako et al. 2001) and NGC 3783 (Kron- 
gold et al. 2003). In all three objects, an 
Fe vil-xm UTA and several Ovil lines are 
the signature of the cooler phase, and the 
Ovill and Fe xvil-xxil are the main fea- 
tures of the hotter phase, with high ioniza- 
tion iron dominating the spectrum below 
15 A. These ionized absorbers are inconsis- 
tent with substantial amounts of material 
with intermediate values of the ionization 
parameter. Furthermore, the ionization 
degree of the LIP is consistent with a cooler 
medium producing significant amounts of 
O VI and other lower ionization species, 
which should also imprint their signature 
in the UV band (see also §5.1). In the 
case of NGC 3783, a third very high ionized 
component has also been detected (Netzer 
et al. 2003). In NGC 985 the presence of 
such component is hinted by the data, al- 
though it cannot be well constrained (see 
Fig. 11). 

The similar characteristics in the ab- 
sorbers of these objects suggest that this 
an intrinsic property related to the struc- 
ture of the nuclear environment. In the 
following sections, a possible scenario re- 
garding the nature of these absorbers is 
discussed. 

5.1. The UV X-ray Connection of 
the Absorbers 

NGC 985 has been observed with HST- 
STIS and FUSE. Arav (2002) modeled 
these data, and showed that, once satura- 
tion is accounted for, the component out- 
flowing at 404 km s" 1 (i.e. component 4) 
produces absorption in the X-ray band. 

In Table 6 we present the O VI, N v, 



C IV, and H I column densities inferred 
from the LIP of our Chandra best fitting 
model. It is evident that this absorbing 
component would produce absorption in 
the UV band. Arav (2002) measured col- 
umn densities of logA^ = 14.2 [cm -2 ], 
and an log Nh { = 14.7 [cm -2 ]. These val- 
ues are smaller than the ones derived from 
our model by factors 29 for N v and 5 for H 
I (see Table 6). Finding consistent column 
densities for the UV and X-ray absorbers 
has been historically difficult as has been 
discussed by Krongold et al. (2003). In 
the case of NGC 985, a few reasons can be 
responsible for the discrepancy: (1) The 
predicted column densities inferred by any 
model in the UV strongly depend on the 
exact shape of the SED assumed in the 
analysis (Kaspi et al. 2001; Steenbrugge 
et al. 2003). Thus, the overestimation in 
the model could be an effect induced by 
our chosen SED (Fig. 3a). In fact, the de- 
rived factors reduce to 15 for N v and 4 for 
H I when considering the SED presented 
in Figure 3b (see Table 6). The reason for 
the decrease in ionic column density in this 
SED is that it has a larger number of pho- 
tons in the extreme UV, the relevant region 
for the abundances of the ions absorbing 
in the UV (note the small decrease in H 
I, since the photon flux assumed in both 
SEDs at 912 A is the same). However, the 
intrinsic SED in this region cannot be con- 
strained due to Galactic absorption. (2) In 
addition to SED effects, the fact that the 
observations are not simultaneous should 
be kept in mind. It is possible that the 
absorber has changed in opacity because 
the gas responds to flux variations (or be- 
cause the absorbing gas is not the same 
at the two different epochs). For instance, 
the ionization parameter for the LIP in the 
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BeppoSAX model is larger than the one 
in the Chandra model. Then, the differ- 
ence in N v column density between Arav's 
measurements and our BeppoSAX best fit 
model is only 13, if the SED presented in 
Fig. 14 is used (this SED is similar in the 
extreme UV to the one presented in Figure 
3a for the Chandra analysis, with a small 
flux in this region). Therefore, the differ- 
ent columns can reflect only the ionization 
state of the gas. (3) Finally, it is possible 
that other UV components (besides com- 
ponent 4) contribute to the absorption ob- 
served in the X-rays. 

Quantifying these differences would re- 
quire simultaneous UV and X-ray obser- 
vations with high resolution and high S/N 
ratio. However, we note that NGC 985, 
along with NGC 5548 (Kaastra et al. 2002; 
Arav et al. 2003; Steenbrugge et al. 2003), 
NGC 3783 (Kaspi et al. 2001; Gabel et al. 
2003; Krongold et al. 2003; Netzer et al. 
2003), and NGC 7469 (Blustin et al. 2003; 
Kriss et al. 2003) are good examples of the 
UV and X-ray absorbing connection. 

5.2. Density and Location of the 
Absorbing Gas 

Our analysis of the LIP indicates that 
this gas may be in photoionization equi- 
librium and responding consistently to 
the changes in the continuum. Using the 
elapsed time between the high state and 
low state observations (1029 days) as an 
upper limit to the recombination time, a 
lower limit on the density of the recom- 
bining gas and an upper limit on the dis- 
tance from the central source can be set 
for this gas. For this calculation, we used 
the photoionization equilibrium temper- 
ature, the flux in the 0.1-10 keV range 



derived from our fit to the Chandra obser- 
vation, the ionization parameter U x , and 
the recombination rate for Ovil (Shull & 
van Steenberg 1982) since it is the dom- 
inant charge state. We estimate a lower 
limit for the electron density of the LIP of 
n e (LIP)>610 cm~ 3 . This translates into 
an upper limit of R(LIP)< 21.6 pc for the 
distance between the central source and 
the ionized gas, and a maximum thickness 
AR(LIP)< 1.3 pc for the gas. 

5.3. Thermal Instabilities and the 
Absorbing Components 

Several models have been proposed to 
explain the nature of the ionized absorbers 
around AGN. The different scenarios in- 
clude a continuous range of ionization pa- 
rameters (Krolik & Kriss 2001), different 
clouds in different regions (including for in- 
stance different UV and X-ray absorbing 
regions), or similar locations for the ab- 
sorbing winds (Krongold et al. 2003). 

The data presented here suggest that a 
continuous range of U does not adequately 
describe the observed spectra (see discus- 
sion in §2.2). The presence of different 
clouds in different regions seems unlikely 
(at least for the objects discussed here), 
not only because of the similar characteris- 
tics of the absorbing components in differ- 
ent AGNs, but also because pressure bal- 
ance is observed (to within the errors) be- 
tween the components (see 5.4). 

An intriguing possibility is that the dif- 
ferent components form because of ther- 
mal instabilities in the thermal photoion- 
ization equilibrium S-curve (Krolik, Mc- 
kee, & Tarter 1982). This curve marks 
the points of thermal equilibrium in the 
T, U/T plane, where T is the equilibrium 
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temperature of the gas, and U/T is a quan- 
tity inversely proportional to the pressure 
of the gas (see Fig. 15). The wiggles 
in the S-curve are due to changing heat- 
ing and cooling rates, as different ions be- 
come dominant. The shape of this curve 
is also affected by the ionizing SED and 
the metallicity of the absorbing gas (e.g. 
Komossa & Fink 1997). The high tem- 
perature branch lies where only Compton 
heating is important. Gas in regions of the 
curve with negative derivative is unstable 
because any isobaric perturbation will be 
amplified, leading to net cooling or heat- 
ing. Then, the different components ob- 
served form simply, when the gas is driven 
to the stable (positive slope) branches of 
the curve, while tending toward equilib- 
rium. 

5.4. Pressure Balance between the 
Absorbing Components 

Krongold et al. (2003) reported that the 
two components in the ionized absorber 
of NGC 3783 are in pressure balance to 
within 7%, which suggested that the ab- 
sorption could arise from two phases of the 
same medium. As we shall point here, this 
also seems to be the situation for the ion- 
ized absorber of NGC 985. 

As can be observed in Figure 15a, the 
HIP and the LIP in the best fit model 
of the NGC 985 Chandra data have very 
different equilibrium temperatures but lie 
close to each other in the U/T axis. As 
mentioned before, with the adopted defini- 
tion of the ionization parameter (U), this 
ratio is inversely proportional to P, the 
gas pressure. Hence, assuming that the 
HIP and the LIP have the same location, 
the gas pressure between them is indistin- 



guishable within 25%, (the pressure of the 
phases is consistent with a single value at 
a level well within 2a). 

We stress again the fact that the absorb- 
ing material does not span all the points 
in the stable branches of the S-curve, but 
rather it is concentrated in separate phases 
(see §2.2). Therefore, although the effect 
of thermal instabilities may be partially re- 
sponsible for the formation of the phases 
(§5.3), this cannot alone explain the ob- 
servations. One possibility is that these 
phases form as the gas tends to go to the 
stable branches of the thermal equilibrium 
curve, but keeping pressure balance be- 
tween them (see below). 

There are two special considerations re- 
garding the validity of the pressure equi- 
librium result that have to be discussed at 
this point: 

(1) The shape of the SED: Kaspi et al. 
(2001) and Steenbrugge et al. (2003) have 
shown that the shape of the SED in the 
far UV (the unobservable region of the 
spectrum) has a negligible effect on the 
charge states absorbing in the X-ray re- 
gion. However, the curve of thermal equi- 
librium does depend on the shape of the 
SED at these wavelengths. Thus, it is pos- 
sible that SEDs with a UV-X ray break dif- 
ferent to the one assumed here would result 
in two phases lying out of pressure balance. 
In other words, our result could be reflect- 
ing specific characteristics of our chosen 
SED. To further investigate this possibil- 
ity, we have repeated our entire analysis 
with the SED presented in Figure 3b (with 
the UV-X ray break at 1 keV instead of 0.1 
keV) . This choice gives the maximum effect 
introduced in our results by SED uncer- 
tainties since the break cannot be at larger 
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energies because then, the observed spec- 
trum would be directly affected. We find 
that pressure equilibrium is found regard- 
less of the assumed shape of the far UV 
SED (see Figure 15b). 

(2) The analysis of the LIP is based mostly 
on the UTA: As discussed in §2.2.1, the 
lack of reliable DRR for low Fe charge 
states may be producing an underestimate 
of the ionization parameter of the LIP. 
Even if this effect is real, it does not mod- 
ify the analysis presented here in any way. 
More accurate DR rates would shift the po- 
sition of the LIP to higher values of both 
pressure and temperature in the S-curve 
(A log 17 < +0.3, see §2.2.1). However, 
the shift is insufficient to move the LIP out 
of pressure equilibrium with the HIP. On 
the contrary, such a change would make 
the pressure of both components even more 
closely matched (Figure 15). 

5.5. Pressure Confinement of the 
Absorbing Components 

A striking result is that pressure equilib- 
rium between the two phases is also found 
in the BeppoSAX observation of NGC 985. 
This can be observed in Figure 16, where 
the S-curve for the SED used in the Bep- 
poSAX analysis is shown, and the positions 
of the HIP and the LIP in this model are 
also marked. Although it is not clear what 
drives the change of ionization in the HIP, 
both phases again align close to a single 
value of logU/T, as would be expected if 
one of the two components confines the 
other one (or as if both components are 
pressure confined by a third phase). 

There are two necessary conditions for 
this to happen: (1) The hydrodynamic 
timescale (t#) in the confined medium has 



to be shorter than the timescale in which 
the gas reacts to the changes of the con- 
tinuum. The sound velocity v s for media 
at temperature ^ 10 6 K is ^ 100 km s _1 , 
and the hydrodynamic timescale can be es- 
timated as t H > AR/v s , where AR is the 
scale size of the medium. For a single gas 
cloud with conditions typical of ionized ab- 
sorbers (10 4 K<T< 10 6 K, and N H > few 
10 21 cm -2 ), the hydrodynamic timescale 
will always be larger than the photoioniza- 
tion timescale, which prevents the gas pres- 
sure from responding to the change in the 
incident continuum quickly enough to re- 
main in pressure equilibrium. This is true 
for any density inferred using the ioniza- 
tion parameters from our model, assuming 
that the gas is located anywhere within 
10 parsecs of the central source. Thus, 
an important conclusion is that the phases 
can only be in true pressure balance if the 
confined medium is composed by a large 
number of small clouds (~ 10 4 clouds), so 
that the hydrodynamic time scale can be 
shorter than the photoionization timescale. 
The total column density of the absorber 
would result from the individual contribu- 
tions of the clouds. 

(2) The photoionization equilibrium 
timescale of the media has to be shorter 
than the variability timescale of the central 
source, in order for the absorbing material 
to be close to photoionization equilibrium. 
This condition is assumed in our analysis. 
Clearly, if the gas is out of photoionization 
equilibrium, then the ionization parame- 
ters and temperatures derived here are not 
meaningful, and the absorbing gas is not 
even constrained to lie on the S-curve. We 
note that photoionization equilibrium of 
the confining medium is not required for 
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pressure confinement (see §5.5.3). 

We will follow the discussion assuming 
that these two conditions are indeed sat- 
isfied (but see 5.6). In such a scenario, 
the confining medium would be respond- 
ing to the changes in the continuum, and 
moving freely on the S- curve. The con- 
fined medium, on the other hand, would 
not be able to move freely in the S-curve, 
rather it would have to follow at all times 
the position of the confining medium in the 
U/T axis. This means that the confined 
medium would have to adapt to the ioniz- 
ing continuum with a constrained value of 
the pressure at all times. Thus, to reach 
the photoionization equilibrium tempera- 
ture, the confined phase would have to ex- 
pand or contract to keep pressure balance 
with the confining medium. Then, the ion- 
ization parameter (and thus, the opacity) 
of the confined phase would be determined 
not only by the change in flux of the cen- 
tral source (as in an isolated system), but 
also by the change in density due to the 
change in volume. 

5.5.1. Pressure Confinement of the HIP 
by the LIP 

This scenario can easily explain the ob- 
served behavior of the HIP between the 
BeppoSAX and Chandra observations: as 
the flux drops, the LIP gets colder and its 
pressure drops. Hence, the HIP would have 
to expand to decrease its pressure. The 
drop in density of the HIP would increase 
its ionization parameter and temperature 
(even though the flux is dropping), until 
photoionization equilibrium is reached. 

However, an inconsistency for this pic- 
ture is that the column density of the LIP 
is smaller than that of the HIP and its den- 



sity is larger (if the phases are at the same 
distance, the difference in U between them 
should be mostly an effect of the differ- 
ent densities). This would require the HIP 
to be more extended than the LIP, and at 
the same time, the HIP would have to be 
made up of several small clouds immersed 
in the LIP. We find it extremely unlikely 
that these two conditions can be satisfied 
simultaneously. 

5.5.2. Pressure Confinement of the LIP 
by the HIP 

A second scenario is that the HIP is con- 
fining the LIP. However, the HIP would 
have to be out of photoionization equi- 
librium, in order to explain the changes 
in ionization parameter. This picture is 
not possible, since as explained before, if 
this phase is not in photoionization equilib- 
rium, then its position on the T-U/T plane 
is unknown. Thus, the alignment on the 
S-curve would have to be considered just a 
bias induced by our analysis (see §5.6). 

5.5.3. Pressure Confinement of the HIP 
and the LIP by a Third Component 

A third scenario is that a more ionized 
component is confining both the HIP and 
the LIP. Being more ionized, this compo- 
nent would require a density lower than 
that of the HIP (which is already 30 times 
smaller than the density of the LIP). In 
this case we have two possibilities. The 
first is that despite the low density, this 
component still has time to adapt to the 
changes in the continuum, and therefore is 
in photoionization equilibrium. The sec- 
ond is that this is not the case. In either 
scenario, the HIP and the LIP would have 
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to follow the position of this hypothesized 
phase in the U/T axis. In the photoioniza- 
tion equilibrium situation the third compo- 
nent has to be located on the curve of ther- 
mal stability. In the non-photoionization 
equilibrium case, the component could be 
located off the curve, as long as the pres- 
sure imposed by this phase on the other 
two is compatible with the region of the S 
curve where multiple phases can survive, 
otherwise the phases would dissolve. 

This scenario is appealing because the 
presence of this phase could easily explain 
the observed change of ionization parame- 
ter of both the LIP and the HIP, between 
the BeppoSAX and Chandra observations, 
as we show below. In addition, the data 
hints the presence of more ionized gas (see 
Fig. 11). In NGC 3783, the presence 
of such absorber has already been estab- 
lished, and found to be in pressure balance 
with the less ionized components (Netzer 
et al. 2003), which further suggests that 
this could be the case in NGC 985. Evi- 
dently, much better data and time evolv- 
ing photoionization models (e.g. Nicastro 
et al. 1999) are needed to explore further 
this scenario. 

If a third confining component is the 
correct picture, then the HIP and the LIP 
have to be close to photoionization equilib- 
rium, but may not respond "as expected" 
to the flux variations because their pres- 
sure is fixed to the pressure of the hottest 
phase. If this hot phase does not change 
its position on the (U/T,T) plane dramat- 
ically in response to continuum changes 
(because it is out of photoionization equi- 
librium or because the gas lies on a steep 
part of the S-curve), then the HIP and/or 
the LIP may appear as if they were not 



following (or not even responding to) the 
continuum. 

5.6. The Absorbing Components 
Are Not in Pressure Balance 

From the last sections, it is clear that a 
scenario of pressure-confining is appealing 
because, due to the regions of thermal sta- 
bility in the S-curve, it can explain the dif- 
ferent absorption components observed in 
the spectra of Seyfert galaxies. However, 
such scenario also imposes strong physical 
conditions on the media that could be dif- 
ficult to satisfy. There is not enough evi- 
dence currently to rule out or confirm such 
a scenario and other possibilities have to 
be considered. 

If the pressure confinement is not real, 
why do the phases align close to a single 
value of U/T? The answer to this ques- 
tion may rely on the shape of the S-curve. 
For several SEDs, the S-curve can be ex- 
tremely steep in a large range of temper- 
atures (between 10 4 to 10 6 K) where pho- 
toionized gas absorbs in the X-ray range 
(Krolik & Kriss 2001). Large gradients of 
temperature (factors 10-50) can then arise 
from small changes of the gas pressure. 
This is clearly illustrated in Figures 15 
and 16. Thus, if several absorbing compo- 
nents were in photoionization equilibrium 
and located in the "vertical region" of the 
S-curve, but were completely disconnected 
from each other, they would appear as if 
they were in pressure balance. While this 
is certainly plausible, we note that there 
are also SEDs that do not show such a pro- 
nounced vertical branch (Komossa & Fink 
1997). In those cases, it is possible to dis- 
entangle true pressure balance from appar- 
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ent pressure balance (see for example the 
S-curve obtained for the SED derived by 
Krongold et al. 2003 in their analysis of 
NGC 3783, their Fig. 16). Thus, discrimi- 
nating true pressure balance from apparent 
pressure balance strongly depends on (1) 
the quality of the X-ray spectrum to con- 
strain better the position of the phases in 
the (T,U/T) plane, and (2) the true shape 
of the SED, which is, however, always un- 
certain. 

Evidently, if the gas is out of photoion- 
ization equilibrium, the curves of thermal 
stability discussed here are not relevant to 
describing the physical state of the gas. 
Then, we would not know if the phases 
are or are not in pressure balance. We 
note that the excellent fits obtained for the 
high quality data of NGC 3783 (Krongold 
et al. 2003, Netzer et al. 2003) suggest 
that the gas is at least close to photoioniza- 
tion equilibrium. Time evolving photoion- 
ization models applied to high quality data 
are the best way to determine if the data 
are in photoionization equilibrium, and to 
further study possible scenarios of pressure 
confinement out of equilibrium. 

Additional complexity is present be- 
cause the X-ray detectors are most sen- 
sitive to absorption features produced by 
gas in exactly the vertical region of the 
S-curve. So, the detection of components 
only in the vertical branch (i.e. in apparent 
pressure equilibrium) could be the reflec- 
tion of an observational bias. To explore 
this possibility we have simulated a series 
of Chandra observations using the same ex- 
posure time and the same best continuum 
model for the NGC 985 observation (see 
Table 2). We have simulated data with 
two different absorption scenarios: 



(1) The first dataset was produced using 
the HIP plus gas less ionized than the LIP 
(i.e. gas located to the left of the LIP 
in the horizontal branch of the S-curve 
with log[/=-0.8 and logN H =21.5 [cm~ 2 ]). 
We modeled these data with PHASE and 
found that the presence of colder gas was 
easily detected and constrained. 

(2) The second dataset was simulated us- 
ing 3 absorption components: the HIP and 
the LIP (Table 3), plus a third cooler ab- 
sorber with relatively large column den- 
sity (log[/=-0.8 and logN H =21.5 [cm~ 2 ]). 
First, we modeled these data including 
only two absorption components. This ap- 
proach gives an acceptable fit, but the ion- 
ization parameter and column density of 
the low ionization gas are only upper lim- 
its, which allows the presence of the ad- 
ditional cold ionized material. Then, we 
fitted the simulated data with three ab- 
sorbing components. With this model the 
values of the column density and ioniza- 
tion parameter of the LIP can be recov- 
ered within a significance of ~ 1.5a, as ob- 
served in Figure 17a, and the parameters 
of the HIP can also be recovered within 
a level ~ la. The third component is 
well constrained in N^ (21.4<N# <21.6), 
but the model is consistent with the ac- 
tual values used to simulate this gas at a 
significance level of only 2.5a, with U 
poorly constrained (to a factor 100). This 
can be observed in Figure 17b, where we 
present the confidence region for this com- 
ponent. In fact, our simulations show that 
the presence of gas 75 times less ionized 
than the LIP (logC/ ~ -2) would still be 
detectable, because the bound-free opacity 
increases dramatically with decreasing U. 
However, the ionization state of such gas 
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would be basically unconstrained (log U < 
-0.5). Thus, even with low S/N ratio data, 
we feel confident that there is no addi- 
tional absorbing material along the hori- 
zontal branch of the S-curve (i.e. gas out of 
pressure balance) in the spectrum of NGC 
985, or it would be detected. Evidently, 
the high quality data of NGC 3783 is less 
susceptible to this kind of bias. 

All these tests indicate that, while the 
pressure-confining scenario might be ap- 
propriate, it requires fragmentation of 
the confined absorber, and thus, the idea 
should be tested even further. Any model 
of the ionized absorbers around AGN has 
to explain several points: (1) Why is there 
not a significant amount of gas on the hor- 
izontal branch of the S-curve, at temper- 
atures lower than the one of the LIP? (2) 
Why does the gas tend to give preference to 
special values of the ionization parameter? 
(3) How do the observed components re- 
late to the stable branches of the S-curve, 
assuming that the gas is indeed in pho- 
toionization equilibrium? (4) What is the 
relation between the UV and the X-ray ab- 
sorbers? In several objects, the observed 
X-ray absorbing components do not oc- 
cupy all the points in the stable branches 
of the S-curve, and the absorption seems to 
be produced from similar charge states. A 
simple model with absorbing components 
located at random cannot explain these 
properties. Rather they seem to point to 
real properties of the structure of active 
nuclei. 

A perhaps more basic question that still 
has to be answered is 'where are the X- 
ray absorbers located?' Are they located 
at large scale distances (parsecs from the 
nucleus; Netzer et al. 2003; Behar et al. 



2003) forming part of the biconical flow 
observed for the narrow emission line re- 
gions (Kinkhabwala et al. 2002; Behar et 
al. 2003)? Or, on the contrary, are they 
closer (at subparsec distances, e.g. Nicas- 
tro et al. 1999; Netzer et al. 2002; Gabel 
et al. 2003; Reeves et al. 2003) and part 
of a disk wind (Elvis 2000)? 

Answers to these questions will help us 
understand what is the genesis of these 
winds and what is their relation to the 
black hole and the accretion disk. 

6. ALTERNATIVE SCENARIOS 
FOR THE IONIZATION MECH- 
ANISM 

The solution found for the ionized ab- 
sorber through photoionization models is 
satisfactory and carries important implica- 
tions for the nuclear environment of AGNs. 
Nevertheless, an alternative to photoion- 
ization by the central source deserves to be 
explored for the HIP, since this phase does 
not respond simply to continuum changes. 

The ring located in NGC 985 at a dis- 
tance of kiloparsecs from the nucleus has 
a large inclination angle, which brings the 
possibility that our line of sight to the 
central source crosses through collisionally 
ionized material in this ring. Such colli- 
sions would be produced by an expanding 
wave triggered by a pole-on merger (see 
for instance Karovska et al. 2002 for Cen- 
taurus A). Hence, there is a chance that 
the material producing the absorption ob- 
served in X-rays is collisionally ionized (the 
HIP is too hot to be photoionized by the 
starburst induced by the merger). If colli- 
sions drive the ionization, the temperature 
of the HIP is at least an order of magni- 
tude hotter than our estimate, and thus 
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the two components cannot be in pressure 
equilibrium. Although this scenario can- 
not be excluded, we notice that the change 
in ionization structure of the HIP observed 
in the 3 year period elapsed between the 
observations, requires an increase in the 
electron temperature by at least a factor 
of 2 in the HIP (the temperature variation 
is significant at a level 2.6a), which sug- 
gests that this gas is probably not ionized 
by collisional mechanisms. Furthermore, if 
the absorption of the HIP arises in the ring, 
then the pressure balance found between 
the phases would have to be interpreted as 
a pure coincidence. We stress again that 
pressure balance is also observed in NGC 
3783, a relatively undisturbed galaxy. 

7. SUMMARY 

A new Chandra observation and an 
archival BeppoSAX observation are used 
to study the nuclear environment of NGC 
985. Absorption features consistent with 
the presence of ionized gas in the nuclear 
environment of this Seyfert 1 galaxy are 
detected in both datasets. We have used 
the code PHASE (Krongold et al. 2003) to 
model in a self-consistent way the ionized 
absorber, and to further study possible 
changes in the opacity of the gas to con- 
tinuum variability. Our main results are 
as follows:. 

(1) The intrinsic continuum of the 
source is well reproduced by a power law 
(T 1.4 — 1.6) and a thermal component 
(kT=0.1 keV). This continuum was atten- 
uated by an equivalent hydrogen column 
density of 3.0 x 10 20 cm -2 to account for 
Galactic absorption. 

(2) The data requires two absorption 
components different by a factor of 29 in 



ionization parameter, and by a factor of 
3 in H column density. The hotter phase 
is clearly identified through broad absorp- 
tion features produced by several blends of 
Fe L-shell absorption, and the cooler com- 
ponent is detected due to the presence of 
a deep Fe M-shell UTA. However, the sig- 
nal to noise ratio of the Chandra data does 
not allow the identification of significant 
(er > 2) narrow (unblended) absorption 
features. This results in a poor determi- 
nation of the outflow velocities of the com- 
ponents (581 ± 206 km s" 1 for the high 
ionization phase and 197 ± 184 km s _1 for 
the low ionization one). 

(3) Despite the limited signal to noise 
ratio of the Chandra data, the ioniza- 
tion state of both components is well con- 
strained thanks to several features of Fe 
L-shell absorption, but mainly because of 
the Fe UTA. The lack of accurate dielec- 
tronic recombination rates of low charge 
states of Fe (Netzer et al. 2003) intro- 
duces an uncertainty by a factor less than 
2 in the ionization parameter of the phases; 
however, our results are robust and highly 
insensitive to this source of error. 

(4) The outflow velocity of both compo- 
nents is consistent with most of the absorp- 
tion components found in the UV spec- 
trum of this object. In particular, both 
components are consistent with the UV 
component with outflow velocity of 404 km 
s _1 , which is hot enough to absorb in the 
X-ray regime (Arav 2002). In addition, 
the low ionization component in the X-rays 
produces significant amounts of O VI, N v, 
and C IV, which strongly suggests that the 
absorption features in both bands are the 
manifestation of a single wind. 

(5) The data set tight limits on the pres- 
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ence of additional cold absorption compo- 
nents, showing that gas with temperature 
< 10 6 K clumps around the two absorbers 
detected. On the other hand, the presence 
of a third component, hotter than the other 
two (T> 10 6 K), and producing Fe K-shell 
absorption, is hinted by the data. How- 
ever, data with better signal to noise ra- 
tio is needed to constrain the properties of 
such component. 

(6) The gas detected in the spectrum of 
NGC 985 has strikingly similar character- 
istics to the gas found in NGC 3783 and 
IRAS 13349+2438: (1) absorption pro- 
duced by cold gas (Fe VII-XIIl), (2) ab- 
sorption by hot gas, dominated by Fe xvil- 
XXII, and (3) no significant absorption con- 
sistent with Fe XIV- XVI. Then a continu- 
ous range of ionization parameters is dis- 
favored (at least in these objects). These 
similar characteristics further suggest that 
this is an intrinsic property of the ab- 
sorbers, probably related to the unstable 
branches of the thermal equilibrium curve. 

(7) The X-ray luminosity drops by a fac- 
tor 2.3 between the BeppoSAX and Chan- 
dra observations (separated by ~ 3 years). 
We detect (model dependent) variability in 
the opacity of the absorbers. While the 
colder component is consistent with a sim- 
ple picture of photoionization equilibrium, 
the ionization state of the hotter compo- 
nent seems to increase while the flux drops. 

(8) The most striking result in our anal- 
ysis is that the two absorbing phases are 
consistent (to within the errors) with pres- 
sure equilibrium. This result is obtained 
during both the Chandra and the Bep- 
poSAX observations. Based on this, we 
speculate that the absorption arises from a 
multi-phase wind. Such a scenario can ex- 



plain the change in opacity of both compo- 
nents during the observations, but requires 
that a third hotter component is pressure- 
confining the two phases. Thus, our anal- 
ysis strongly points to a 3-phase medium, 
just like the wind suggested in NGC 3783. 

(9) An important implication of our 
analysis is that, if the pressure-confining 
scenario is real, then fragmentation of the 
confined phases in a large number of clouds 
is required. 

The analysis presented in this paper 
provides further evidence that warm ab- 
sorbers in AGN are simple systems that 
can be relatively well understood. The 
high resolution spectroscopic capabilities 
of the new generation of X-ray observa- 
tories have revealed warm absorbers as 
strong winds, now recognized as an im- 
portant feature of AGN structure. Here, 
we have used high and low resolution data 
to constrain and characterize the physical 
conditions and structure of these winds. 
Measuring the response in the opacity to 
variability of the ionizing flux, we suggest 
that the nature of the ionized absorber is 
that of a multi-phase wind, with one phase 
confining the other two. The implications 
of this result are extremely important to 
our understanding of quasar energy bud- 
gets, as well as to our understanding of 
the interaction between these winds and 
the interstellar medium of the host galaxy. 
However, our analysis is based on a com- 
parison between high and low resolution 
data. Only a direct comparison between 
high resolution data will allow further test- 
ing of the scenario presented here, leading 
to a deeper understanding of the nature of 
AGN and quasars. 
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Table 1 
Observation Log of NGC 985 



Observatory 


Sequence Number 


UT Start 


Exp Time (ks) a 


BeppoSAX 


50862001 


1999 AUG 20 06:38 


95.0 


Chandra 


700449 


2002 JUN 24 19:33 


77.7 



Sum of good time intervals. 
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Table 2 
Continuum Parameters 



Parameter 


Chandra Data 


BeppoSAX Data 




0.7-8 keV a 


3-8 keV 


0.1-110 keV a 


r b 


1.60±0.03 


1.62±0.02 


1.40±.04 


Norm c 


1.22±0.04 


1.22±0.04 


2.47±0.04 


BB kT d 


0.10±0.01 




0.10±0.02 


Norm 6 


3.7±0.8 




8.8±1.2 


E c (keV)' 






67 ±24 


Flux(.l-2keV) s 


0.63±0.12 




1.34±0.08 


Flux(2-10keV) s 


0.54±0.09 




1.49±0.12 


R h 






0.0+0.35 


X 2 /d.o.f. 


179.0/174 


46.7/49 


116.4/110 



a Continuum fits carried out with the self-consistent inclusion 
of an ionized absorber. See Tables 3 and 4. 

b Power Law Photon Index 

Tower law normalization in 10 -3 photons keV -1 cm~ 2 s" 1 
at 1 keV 

d Blackbody Temperature in keV 

c Blackbody Normalization in 10~ 5 L 3 g/Df , where L 39 is the 
source luminosity in units of 10 39 erg s _1 and D w is the distance 
to the source in units of 10 kpc 

f High Energy Cutoff as inferred from the SAX data 

g In units of 10 -11 ergs cm~ 2 s _1 

h Reflection component factor (R=Afi/27r, with Af2 the solid 
angle subtended by the reflector). 
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Table 3 

Two Phase Absorber Parameters for Chandra Data 



Parameter 



High-Ionization Low-Ionization 



HIP LIP 



Log U a 

Log Ntf (cm- 2 )" 
V Turb (km s- 1 ) 6 
Vout (km s- 1 )" 
[Log T (K)] c 
Log U/T (ocl/P d ) 
Log U« 



1.34±0.10 
21.81±.22 



21.39±0.19 
350 
197±184 
4.47±0.04 




350 



581±206 
5.79±0.11 
-4.45 ±0.21 
-1.09±0.10 



-2.55±0.09 




a Free parameters of the model. 
b Fixed parameter. 

c Derived from the column density and ionization pa- 
rameter, assuming photoionization equilibrium. 

d The gas pressure Pocn e T. Assuming that both phases 
lie at the same distance from the central source n e ocl/U, 
and PocT/U. 



ionization parameter defined in the 0.1-10 keV range 
(Netzer 1996). 
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Table 4 

Two Phase Absorber Parameters for BeppoSAX Data 



Parameter 


High-Ionization 


Low-Ionization 




HIP 


LIP 



Log U a 0.83±0.13 -0.22+°;3° 

Log Ntf (cm- 2 ) 6 21.81 21.39 

V Turb (km s- 1 ) 6 350 350 

Vout (km s- 1 ) 6 581 197 

[Log T (K)] c 5.48±0.12 4.56±0.05 

Log U/T (ocl/P d ) -4.65 ±0.25 -4.78±g"?g 

Log W x -1.25±0.13 -2.30±0.10 



a Free parameter of the model. 

b Fixed parameters, assuming N#, FWHM, and Vout 
values from the best fitting Chandra model. 

c Derived from the column density and ionization pa- 
rameter, assuming photoionization equilibrium. 

d The gas pressure Pocn e T. Assuming that both phases 
lie at the same distance from the central source n e ocl/U, 
and PocT/U. 

ionization parameter defined in the 0.1-10 keV range 
(Netzer 1996). 
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Table 5 

Ionic Column Densities' 1 predicted by the Models 



HIP LIP 



Ion 


Chandra 


BeppoSAX 


Ion 


Chandra 


BeppoSAX 


P \n 
U VI 


1 P K1 
ID. Do 


10. oo 


P \r 
Kj V 


1 i .61 


I ( .Ul 


U VII 


lo.o4 


1 8 IP 
lo.oO 


P \n 
U VI 


17 PA 
I ( .04 


1 ( .Oo 


IN VII 


1 P A P 
10.40 


1 p an 

10. oU 


IN VI 


1 7 

I ( .Uo 


IRQ/I 
10. »4 


N vin 


17.76 


17.73 


N vn 


16.95 


16.97 


vii 


16.18 


16.94 


O vn 


18.03 


18.07 


vin 


17.75 


18.09 


O vin 


17.48 


17.69 


ix 


18.61 


18.54 


Ne VI 


17.16 


16.97 


Ne ix 


16.52 


17.09 


Ne vn 


16.70 


16.77 


Ne x 


17.38 


17.54 


Ne vin 


16.44 


16.75 


Ne XI 


17.66 


17.45 


Ne ix 


16.20 


16.78 


Mg xi 


16.55 


16.99 


Mg vin 


16.51 


16.21 


Mg xii 


17.10 


17.02 


Mg ix 


16.54 


16.57 


Mg xiii 


16.91 


16.56 


Mg x 


16.10 


16.45 


Si xiii 


16.96 


17.09 


Si vin 


16.29 


15.76 


Si xiv 


16.98 


16.79 


Si ix 


16.52 


16.34 


Si xv 


16.47 


16.03 


Si x 


16.37 


16.53 


Fe xvn 


16.16 


16.60 


Fe ix 


16.01 


15.28 


Fe xvm 


16.69 


16.83 


Fe x 


16.24 


15.77 


Fe xix 


16.81 


16.70 


Fe XI 


16.34 


16.11 


Fe xx 


16.72 


16.36 


Fe XII 


16.19 


16.33 


Fe xxi 


16.28 


15.68 


Fe xiii 


15.92 


16.28 



log of column density in cm 
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Table 6 

Column densities' 1 of low ionization LIP ions from the Chandra best fit 

MODEL. 



Ion 


SEDl 6 


SED2 C 


VI 


17.23 


17.13 


N v 


15.67 


15.39 


C iv 


15.82 


15.23 


H I 


15.37 


15.33 



a log of column density 
in cm~ 2 

b SED presented in Fig. 

3a 

a SED presented in Fig. 

3b 
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Fig. 1. — Fluxed, combined (MEG+HEG) first-order spectrum of NGC 985. (a) Continuum fitted to the 
spectrum. The three solid lines stand for (from upper to lower): intrinsic continuum of the source, contin- 
uum attenuated by Galactic absorption, and continuum further attenuated by the bound-free photoelectric 
absorption edges. The dotted line represents the predicted power law without the contribution from the 
blackbody component, (b) Two phase absorption model plotted for comparison. Typical errors (per bin) are 
marked with green solid ranges. [See the electronic edition of the Journal for a color version of this figure] 
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Fig. 2. — Fit residuals for continuum determinations over Chandra HEG+MEG data of 
NGC 985. All models are attenuated by cold gas with logN# = 20.5 [cm~ 2 ] due to Galactic 
absorption, (a) Model including a continuum power law fitted over the entire spectral range 
(0.6-8keV). (b) Power law fitted only in the 3-8 keV range, (c) Model in Panel (b) plus 
a soft excess parameterized with a blackbody component. The data are binned to > 100 
photons in each channel to make evident the overall features. [See the electronic edition of 
the Journal for a color version of this figure] 
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Fig. 3. — Spectral Energy Distribution used to model the ionized absorber in the Chandra- 
HETG data. Filled circles (blue) mark the observed SED of NGC 985 (obtained from NED). 
The plot shows the slope adopted in each energy range. This slope relates to the photon 
index as T=s+2. (a) Model with low energy far UV cutoff at 0.1 keV and including the 
contribution of a blackbody component, (b) Model with high energy far UV cutoff at 1 keV. 
[See the electronic edition of the Journal for a color version of this figure] 
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Fig. 4. — Two phase absorber model plotted against the first-order combined MEG+HEG 
spectrum of NGC 985. Absorption lines predicted are marked along bars at top for each 
element (red). The continuum level (including edge continuum absorption) is overplotted 
for comparison (dotted green line). The spectrum is presented in the rest frame system of 
the absorbing gas. Below 10 A the data are presented in bins of size 0.05 A, above 10 A in 
bins of size 0.1 A. [See the electronic edition of the Journal for a color version of this figure] 
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Fig. 5. — Upper panel: MEG+HEG spectrum of NGC 985 grouped in bins of size 0.15 A (b inning factor 
of 30). The solid line (green) marks the continuum level, including photoelectric bound- free absorption from 
our model. The final model including resonant absorption is marked with the dotted line (red) . The feature 
in the continuum (solid green line) between 13.3 and 13.8 A is produced by a chip gap in the detector. 
Lower panel: Residuals (data minus continuum green solid line) in units of standard deviations showing 
the presence of deep Fe L-shell resonant absorption (significance « 3a per bin) and Fe M-shell inner shell 
resonant absorption (significance w 2a per bin). [See the electronic edition of the Journal for a color version 
of this figure] 
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Fig. 6. — Upper Panel: Resonant bound-bound absorption produced by the low ionization 
phase (LIP) of our model, plotted against the MEG+HEG spectrum of NGC 985 grouped 
in bins of size 0.15 A (b inning factor of 30). The dotted line (green) marks the continuum 
level, including bound-free photoelectric absorption from our model. The main feature of 
this component is the Fe UTA. Lower panel: As upper panel, but showing the contribution 
from the high ionization phase (HIP). Below 15 A the spectrum is dominated by Fe L-shell 
absorption. [See the electronic edition of the Journal for a color version of this figure] 
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Fig. 7. — Theoretical transmission spectrum of the two absorption components at 0.001 A 
resolution. The low ionization component produces most of the bound-free photoelectric 
absorption and the UTA near 16 A. The high ionization component produces the Fe L-shell 
absorption below 15 A. [See the electronic edition of the Journal for a color version of this 
figure] 
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Fig. 8. — Outflow velocities for the HIP and the LIP, and for the 6 absorption components 
observed in the UV. The HIP is consistent within 2a with UV components 1 to 5, and the 
LIP with components 3 to 6. Both the HIP and the LIP are consistent with component 4, 
which was expected to absorb in the X-ray region (Arav 2002). 
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b) More Ionized LIP (logU = .4) 
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d) Third Absorber (logU = .6: logN H = 21.5) 
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Fig. 9. — Different values for the ionization parameters of our two phase absorber showing inconsistencies 
with the observed UTA. The left panels present the data in the range 11-14.5 A in bins of size 0.15 A (b inning 
factor of 30), the right panels the data in the range 14.5-19 A in bins of size 0.25 A (b inning factor of 50). 
The green dotted line shows the best fit model for comparison. Panel (a): Less Ionized HIP, logU=0.8. Note 
that the absorption from Fe L-shell is underestimated. Panel (b): More ionized LIP, the ionization parameter 
of the low ionization component has been set to logU=0.4. Panel (c): Less Ionized LIP, logU=-0.8. Panel 
(d): A third absorber has been included in our model with logU=0.6 and logN H=21.5 [cm -2 ]. Including 
a new component overpredicts the absorption by the UTA, which further suggests that the absorbing gas 
tends to clump. [See the electronic edition of the Journal for a color version of this figure] 
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Fig. 10. — Confidence regions constraining the presence of a third absorber in the Chandra 
spectrum of NGC 985. The solid lines correspond to the la, 2a, and 3a levels. Panel (a): A 
low ionization absorber with low column density (logU=-1.0 and logNH=20.8 [cm -2 ]). Panel 
(b): A high ionization component with high column (logU=2.1 and logNH=21.8 [cm -2 ]). 
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Fig. 11. — Residuals to our best fit model over Chandra data in the ranges 1.1-2 A and 10.2- 
11 A. The presence of absorption by material with significant amounts of Fe xxiv-xxvi is 
suggested in the figure, but better S/N ratio is required to constrain the physical properties 
of such component. [See the electronic edition of the Journal for a color version of this 
figure] 
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Fig. 12. — Fit residuals for a model including only a continuum power law constrained 
the 3-8 keV range and further attenuated by Galactic absorption, over BeppoSAX data 
NGC 985. [See the electronic edition of the Journal for a color version of this figure] 
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Fig. 13. — BeppoSAX spectra of NGC 985, plotted against an ionized absorber model 
for comparison, (a) Fluxed spectrum. The blue dashed line shows the continuum further 
attenuated by Galactic absorption, (b) Empirical Spectrum. [See the electronic edition of 
the Journal for a color version of this figure] 
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Fig. 14. — Spectral Energy Distribution used to model the ionized absorber observed in the 
BeppoSAX spectrum. Filled circles (blue) mark the observed SED of NGC 985 (obtained 
from NED). The plot shows the slope adopted in each energy range. This slope relates to 
the photon index as T=s+2. The model assumes a low energy far UV cutoff at 0.1 keV and 
includes the contribution of a blackbody component. [See the electronic edition of the Journal 
for a color version of this figure] 
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Fig. 15. — Curve of thermal stability for the two different SEDs used in the analysis of 
Chandra data, as described in the text, (a) Far UV cutoff at 0.1 keV plus blackbody (Fig. 
3a). (b) Far UV cutoff at 1 keV (Fig. 3b). The positive error bar of the LIP represents 
the maximum effect in U induced by the lack of low temperature DRR. The HIP and the 
LIP are consistent with pressure equilibrium. [See the electronic edition of the Journal for 
a color version of this figure] 
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Fig. 16. — Curve of thermal stability for the SED used in the BeppoSAX analysis, as 
described in the text. The positive error bar of the LIP represents the maximum effect in U 
induced by the lack of low temperature DRR. Pressure balance between the LIP and the HIP 
appears to be present in the BeppoSAX data, even though the ionization degree of the gas 
is different than the one obtained for Chandra data (see text for details). [See the electronic 
edition of the Journal for a color version of this figure] 



46 



21.5 



21.45 



21.4 



21.35 



21.3 



(a) LIP 




21.4 



(b) 3rd Component 




-0.5 



0.5 



-1.5 



-0.5 



logU 



logU 



Fig. 17. — U vs. N# confidence regions for a model fitted to simulated data including 3 
absorption components: the HIP and LIP from the Chandra best fit model plus a third colder 
absorption component with logU=-.8 and logN#=21.5 [cm -2 ]. The solid lines correspond to 
the la, 2a, and 3a levels. In both panels the crossed-circle point marks the value used to 
produce the simulated data and the cross point marks the value obtained in the model. Panel 
(a): Confidence region for the LIP. Panel (b): Confidence region for the third component. 
Clearly, a third component would be detectable in the data, however, the ionization state of 
such a component would be difficult to constrain. 
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